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Symmetry-breaking instability of a laterally uniform granular cluster (strip state) in a prototypical 
driven granular gas is investigated. The system consists of smooth hard disks in a two-dimensional 
box, colliding inelastically with each other and driven, at zero gravity, by a "thermal" wall. The 
limit of nearly elastic particle collisions is considered, and granular hydrodynamics with the Jenkins- 
Richman constitutive relations is employed. The hydrodynamic problem is completely described 
by two scaled parameters and the aspect ratio of the box. Marginal stability analysis predicts a 
spontaneous symmetry breaking instability of the strip state, similar to that predicted recently for 
a different set of constitutive relations. If the system is big enough, the marginal stability curve 
becomes independent of the details of the boundary condition at the driving wall. In this regime, the 
density perturbation is exponentially localized at the elastic wall opposite to the thermal wall. The 
short- and long-wavelength asymptotics of the marginal stability curves are obtained analytically in 
the dilute limit. The physics of the symmetry-breaking instability is discussed. 

PACS numbers: 45.70.Qj 



I. INTRODUCTION 

Granular materials play an important role in indus- 
trial applications, geophysics and astrophysics. They are 
also of a great general interest to physicists, as each of the 
"phases" of granular materials: solid, liquid and gas have 
unusual properties that distinguish them from their clas- 
sic (atomic) counterparts |l], ^|. We will consider gran- 
ular gas (or rapid granular flow) and focus on a variant 
of clustering instability, a striking tendency of granular 
gases to form dense clusters. Clustering instability was 
first discussed in the context of a freely " cooling" granu- 
lar gas j|, ||, |(J. Following these works, related clustering 
phenomena were investigated in driven granular gases as 
well, both in experiments |^] and in particle Simula- 
tions 

Granular clustering results from energy losses by in- 
elastic collisions, and it is a manifestation of thermal 
condensation instability, also encountered in other fields, 
for example in gases and plasmas that cool by their own 
radiation fl^| . Since the discovery of the clustering in- 
stability, the validity of granular hydrodynamics (see Ref. 
|jl3f for a review) has been under scrutiny ||] . In contrast 
to the clustering in a freely "cooling" granular gas, where 
one deals with a complex time-dependent process, steady 
states are achievable in driven granular systems. One of 
the simplest settings of this type is driving the granulate 
by a side wall at a zero gravity. Therefore, an ensemble of 
inelastically colliding hard spheres, confined in a box and 
driven by one or two "thermal" walls has served as a pro- 
totypical driven granular system [jj], §, H> [0j> ■ 
Steady states of this system have served as test beds 
for granular hydrodynamics and its violations. The first 
analysis of this system in the physical literature was per- 
formed, in one dimension, by Kadanoff et al. 0. The 
nearly elastic particles were constrained to move on a 
straight line with energy input from the boundaries. Par- 



ticle simulations |9| showed that, for typical initial con- 
ditions, the system evolves to a state where the particles 
are separated into two groups. Almost all particles form 
a cluster in a small region of space, where they move with 
very small velocities, while a very few remaining particles 
move with high velocities. Clearly, this steady state can- 
not be described by granular hydrodynamics (actually, 
by any coarse-grained theory). Therefore, the results of 
Kadanoff et al. Q brought into question the validity of 
granular hydrodynamics in general. 

This question was addressed in two subsequent the- 
oretical works |l(], [ll| that dealt with similar systems 
in two dimensions. Esipov and Poschel [Tl|| investi- 
gated an ensemble of nearly elastic hard disks in a cir- 
cular box with the circumference serving as thermal 
wall. Grossman et al. [ [To[ considered a rectangular 
box, one side of which served as thermal wall. Parti- 
cle simulations [llj showed granular clusters: dense 
and "cold" regions of granulate developing away from the 
thermal wall. In terms of thecoarse-grained particle den- 
sity, these steady-state clusters had simple shapes: az- 
imuthally uniform (circular state) |Tl| ] and laterally uni- 
form (strip state) |t0|. Grossman et al. also showed 
that, for nearly elastic collisions, the strip state is de- 
scribable by a steady-state solution of granular hydro- 
dynamic equations. The empiric constitutive relations 
suggested by Grossman et al. used simple interpolations 
between the low-density limit, where the constitutive re- 
lations are derivable systematically Jl^] , and high-density 
limit where, close to the dense packing, free volume ar- 
guments can be used. 

The results of Refs. [[t0| [llj showed that the anomaly 
observed in the one-dimensional setting ^ does not per- 
sist in higher dimensions. Clustered states qualitatively 
similar to those of Grossman et al. were observed in ex- 
periment of Kudrolli et al. who investigated a system 
of spherical particles in a box, rolling on a smooth sur- 
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face and driven by a rapidly vibrating side wall. The 
number of particles served as the control parameter in 
Ref. 0. A dense cluster of the strip type was observed 
when the number of particles was big enough, in much 
the same way as in Ref. |To| . The basic physics of the 
strip state is simple and can be explained by the fol- 
lowing hydrodynamic argument. Because of the inelastic 
collisions the particle random motion slows down (that 
is, the granular temperature decreases) with the increase 
of the distance from the driving wall. To maintain the 
momentum balance, the granular density should increase 
with this distance. When the total number of particles is 
big enough (the rest of parameters being the same), the 
density contrast becomes large, and the enhanced den- 
sity region away from the driving wall is observed as the 
strip state. 

The prototypical system exhibits many interesting 
phenomena of non-hydrodynamic nature. These include 
inelastic collapse pL , possible lack of scale separation 
[10 1, non-Gaussianity in the particle velocity distribution 
[10, |R3), normal stress difference and pressure nonunifor- 
For nearly elastic collisions, however, 
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etc. 



granular hydrodynamics was shown to yield an accurate 
quantitative description in the dilute limit |l4| ] , and rea- 
sonably accurate description for moderate and high gran- 
ular densities Of course, the nearly elastic limit is 
quite restrictive for most of realistic granular flows. Still, 
this limit is conceptually important just because gran- 
ular hydrodynamics can be used there. Granular hy- 
drodynamics has a great predictive power and helps to 
identify important collective phenomena (shear flows and 
vortices, shocks, different modes of clustering flows etc.) 
that are difficult, if not impossible, to identify and predict 
in the language of individual particles. Once identified, 
these phenomena can then be investigated in experiment 
and simulations in more general (not necessarily hydro- 
dynamic) formulations. 

Therefore, granular hydrodynamics provide a leading- 
order approach to a big ensemble of nearly-elastically col- 
liding grains. This approach has been taken recently by 
Livne et al. [ p"7| who employed granular hydrodynamics 
for a stability analysis of the strip state described above. 
The analysis revealed a spontaneous symmetry-breaking 
instability of the strip state with respect to perturbations 
along the strip. Well within the instability region, the 
clustered states with broken symmetry, found by a nu- 
merical solution of the steady state hydrodynamic equa- 
tion, are strongly localized in the lateral direction: most 
of the particles are located in dense "islands" jlTj . These 
results indicate that the prototypical system can show a 
non-trivial behavior even in the leading-order, hydrody- 
namic limit. Indeed, this systems can be put into the list 
of pattern- forming systems far from equilibrium Jl9| ] . 

The present work focuses on a more detailed stability 
analysis of this system. Our first objective is to check 
to what extent the symmetry-breaking instability pre- 



dicted in Ref. [[L7| is sensitive to the constitutive rela- 
tions. Livne et al. |l7j employed the empiric relations 
suggested by Grossman et al. jl(|. Here we shall use 
the better known Jenkins- Richman (JR) relations p0[ . 
While the relations of Grossman et al. are more accurate 
for high densities (even including those close to the dense 
packing limit), the JR relations should work better at low 
and intermediate densities. We shall see, however, that 
the marginal stability curves, obtained with these two 
sets of relations, are not much different from each other. 
This implies that the symmetry-breaking instability is 
robust. Our second objective is to get more insight into 
the marginal stability problem and, where possible, to 
obtain analytic results. We shall show that the marginal 
problem is equivalent to an eigenvalue problem of quan- 
tum mechanics. An important finding here is a univer- 
sal behavior of the marginal stability curves in the limit 
when the density perturbations are strongly localized at 
the elastic wall opposite to the thermal wall. In the di- 
lute limit, we obtain analytically the short- and long- 
wavelength asymptotics of the marginal stability curves 
and density eigenfunctions. We also give the physical in- 
terpretation to the symmetry-breaking instability and to 
the density borders of the instability region. 

The rest of the paper is organized as follows. In Sec. 
2 we formulate the model and briefly describe the strip 
state: the simplest steady state of the prototypical sys- 
tem. Section 3 presents marginal stability analysis of 
the strip state and compares the results obtained for two 
different sets of constitutive relations. More results on 
marginal stability, including some analytic results in the 
dilute limit, are presented in Sec. 4. Section 5 includes 
a discussion and summary. 



II. PROTOTYPICAL SYSTEM AND STRIP 
STATE 



The prototypical driven granular system in two dimen- 
sions include inelastically colliding hard disks of diame- 
ter d and mass m = 1, moving in a box with dimensions 
L x x L y . The gravity force is zero. Collisions of disks 
with the walls x = 0, y — and y — L y are assumed 
elastic. The wall x = L x is "thermal" wall: upon colli- 
sion a particle is assigned a random velocity taken from 
a Gaussian ensemble with temperature To- Energy input 
at the thermal wall balances the energy dissipation due to 
inter-particle collisions, so the system can reach a steady 
state. We shall parameterize the inelasticity of the parti- 
cle collisions by a constant normal restitution coefficient 
r and work in the nearly elastic limit: 1 — r 2 <C 1. In 
this limit, the Navier-Stokes granular hydrodynamics is 
expected to be sufficiently accurate in a system with a 
big number of particles and small Knudsen number. The 
possible steady states of the system are described by the 
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steady state versions of the momentum and energy bal- 
ance equations: 



p = const and V • (kVT) = I . 



(1) 



where p is the granular pressure, T is the granular tem- 
perature, k is the thermal conductivity and / is the 
rate of energy losses by collisions. We assume that the 
number density n is not too big: n/n c < 0.5, where 
n c = 2/{\fid 2 ) is the (hexagonal) dense packing density. 
This assumption enables us to employ the constitutive 
relations derived by Jenkins and Richman [^o| . For the 
steady state problem, the required constitutive relations 
include the equation of state p = p{n, T) and relations 
for k and / in terms of n and T. 

Let us introduce scaled coordinates: r/L x — ► r. In the 
new coordinates the box dimensions are 1 x A, where A = 
L y /L x is the aspect ratio. Introducing the normalized 
inverse density z(x,y) = n c /n(x,y), one can rewrite the 
energy balance equation in Eq. (fy) in terms of z(x,y): 



V-(F(z)Vz)=r,Q(z), 



(2) 



where F(z) = A(z) B(z), 
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and rj = (27r/3)(l — r)(L x /d) 2 . Notice that, for an arbi- 
trary small but finite inelasticity 1 — r, the dimensionless 
parameter 7/ can be made arbitrary large, if the system 
size L x is large enough. The parameter n differs from the 
parameter C used by Livne et al. [T^ ] only by a numer- 
ical factor of order unity. Of most interest are regimes 
where r\ 3> 1, see below. 

The boundary conditions for Eq. (g) are determined 
by the properties of the particle-wall interactions. At the 
elastic walls x — 0, y = and y = A we should prescribe 
a zero normal component of the heat flux. In terms of the 
inverse density z we have V„z = at these three walls. 
Here index n denotes the gradient component normal to 
the wall. The constant temperature at the "thermal" 
wall x — 1 yields the condition dz(x = l,y)/dy = 0. To 
make the formulation of the problem complete, one more 
condition is needed. In experiment or particle simula- 
tions, the number of particles N is fixed. This yields a 
normalization condition: 
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where / = (n) jn c is the area fraction of the grains and 
(n) = N/(L x L y ) is the average number density of the 
grains. 

Equations (|2])-(fi]) and the boundary conditions make 
a complete set. One can see that the governing param- 
eters of this system are the scaled parameters rj, / and 
A. If the system is infinite in the y-direction, only two 
governing parameters: n and / remain. Notice that the 
steady-state density distributions are independent of the 
wall temperature To, in contrast to the similar problem 
with gravity, where the gravity acceleration, combined 
with To and the (finite) system size in the direction of 
gravity, would make an additional governing parameter. 

The laterally uniform steady state (strip state) corre- 
sponds to the one-dimensional (y-independent) solution 
z = Z(x). It is described by the equations 



(FZ'Y = r)Q, Z'\ x=0 = 0, 

and / Z~ (x)dx = f, 
Jo 



(5) 



where primes stand for the x derivatives. For the strip 
state, the boundary condition at the wall x = 1 drops 
out. This implies, in particular, that the density profile 
of the strip state is independent of the exact nature of 
the driving wall (thermal or vibrating wall) 0]. This 
degeneracy of the strip state is caused by the character 
of particle interaction: the hard-core potential does not 
introduce any characteristic energy Q. Notice that, in- 
stead of prescribing the grain area fraction /, one can 
prescribe the inverse density Z = Zq at x = 0. This 
condition, combined with the no-flux condition at x = 
defines a Cauchy problem for Z(x). Solving the Cauchy 
problem, one can then compute, from the last equation 
in Eq. (||), the respective value of /. At fixed n, there is a 
one-to-one correspondence between Zq and /. Therefore, 
an alternative parameterization of the strip state is given 
by the scaled numbers n and Zq. We shall see below that 
the same property keeps (and can be conveniently used) 
in the marginal stability problem. Figure [l] shows a typ- 
ical example of the scaled density profile n(x)/n c of the 
strip state obtained by solving Eqs. (||) numerically. 



III. INSTABILITY OF THE STRIP STATE: 
MARGINAL STABILITY CURVES 

In general, the strip state is only one of the possible 
solutions of Eq. (^) . Because of its nonlinearity, Eq. (^) 
may have additional solutions satisfying the same bound- 
ary conditions. When exist, these additional solutions 
are truly two-dimensional: the translational symmetry 
along y is broken. An important class of these solutions 
bifurcate supercritically from the strip state jl7) . There- 
fore, close to the bifurcation point, these solutions can 
be found by linearizing Eq. (0) around the strip state. 
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FIG. 1: An example of the scaled density profile n(x)/n c of 
the strip state for the JR constitutive relations. The governing 
parameters are n = 10 4 and / = 0.0342. 



A similar analysis was performed in Ref. |17[ for the 
constitutive relations of Grossman et al. |10( |. In the 
framework of a time- dependent hydrodynamic formula- 
tion, this analysis corresponds to marginal stability anal- 
ysis of the strip state with respect to small perturbations 
along the strip [jlTj . 

Substituting z(x,y) = Z(x) + ipk{x) cos ky and lin- 
earizing Eq. (j^) with respect to the small correction 
ipk{x) cos ky, we obtain: 
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Here <f>(x) — F(x) ipk(x), functions F and Q are evaluated 
at z = Z{x), and subscript Z means the z derivative 
evaluated at z = Z(x). The boundary conditions are 



b'(x = 0) = and 4>{x = 1) = . 



(7) 



Equation (|^) coincides with the Schrodinger equation for 
an even wave function 4>{x) of a particle in the potential 
well 



U(x) = i ? ^ W < 1 ' 

' L oo otherwise . 



(8) 



The quantity — k 2 serves as the energy eigenvalue. The 
energy levels in the potential (||) are always discrete, and 
there is an infinite number of them. However, as the 
wave number k should be real, only negative or zero en- 
ergy levels are admissible. At fixed values of r\ and /, 
the potential (||) admits at most one such energy level. 
The absence of negative energy levels implies that, in the 
vicinity of the strip state, there are no steady states dif- 
ferent from it. The presence of a negative energy level 
corresponds to a "weakly two-dimensional" steady state, 
bifurcating from the strip state. We shall exploit the 
quantum-mechanical analogy more fully in Sec. 4. Here 
we report some numerical results. Figure shows the 



marginal stability curves: the curves k — k(f) at differ- 
ent values of t], computed numerically. In these compu- 
tations, the parameter r\ was taken large enough. The 
strip state is unstable below the respective curve and 
stable above the curve. Notice that, at fixed rj, the in- 
stability is possible only within a finite interval of /: 
fi(v) < / < fii 7 ])- The same property was reported 
in Ref. |f7j for another set of constitutive relations. We 
shall give a physical explanation to this finding in Sec. 
4. Notice (see also Ref. Q) that, at large rj, the high- 
density stability border ]i is quite small. The curves in 
Fig. H are actually plotted in scaled coordinates: krf~ x l 2 
versus fr/ 1 / 2 . It can be seen that, in the scaled coor- 
dinates, all the curves exit from the same point of the 
horizontal axis frj 1 ^ 2 . In addition, the maxima of all 
the curves are equal. These observations will be also ex- 
plained in Sec. 4. 
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FIG. 2: Marginal stability curves for different values of n, 
plotted in scaled coordinates: kn~ ' versus fn 1/2 . For a 
fixed n the strip state is stable above the respective curve and 
unstable below the curve. The values of n are: 10 4 (solid line), 
2.5 • 10 4 (dashed line), 5 • 10 4 (dotted line) and 10 5 (dash-dot 
line) . 

If the system is infinite in the lateral direction, A = oo, 
while rj and / are fixed, a whole continuum spectrum of 
wave numbers between k = and k = k(rj, f) is unstable. 
Both in experiment, and in numerical simulations A is 
finite. In this case k becomes discrete because of the 
boundary conditions: k = mir/A, where m = 1,2,.... 
For each m we can find the critical value of the aspect 
ratio A (let us call it A m ) such that for A > A m the 
strip state looses stability with respect to the m-modc. 
Obviously, A m = mAi, and Ai is the lowest critical 
value for the symmetry-breaking instability. Figure ^ 
shows Ai as a function of / for different values of 77. 
For fixed 77 and /, the strip state is unstable above the 
curve. It is seen from Fig. [3] that, in order to observe 
the symmetry-breaking instability, one does not need to 
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work with very large aspect ratios A: it is sufficient if the 
system is big enough, so that parameter r\ is sufficiently 
large. 
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FIG. 3: The critical aspect ratio Ai for the symmetry- 
breaking instability as a function of / for different values of 
rj. For a fixed r\, the strip state is stable below the respec- 
tive curve and unstable above the curve. The parameters are: 
771 = 10 4 (solid line), 772 = 2.5 ■ 10 4 (dashed line), 773 = 5 ■ 10 4 
(dotted line) and 774 = 10 5 (dash-dot line). 



To what extent is the symmetry-breaking instability 
sensitive to the precise form of the constitutive relations? 
We compared the marginal stability curves A = Ai(/) 
for different values of 77 with the respective curves ]Tij ] 
found for the constitutive relations of Grossman et al. 
pXj| . A typical example of this comparison is shown in 
Fig. ||. One can see that, qualitatively, the results are 
the same: the both curves describe a symmetry-breaking 
instability at a critical value of the aspect ratio that 
depends on the area fraction. In both cases, there are 
sharp low- and high-density borders of instability region. 
Therefore, we can conclude that the instability is robust 
and does not require a very special form of the consti- 
tutive relations. On the other hand, there is a notice- 
able (about 15%) difference in the exact positions of the 
marginal stability curves, so the instability provides a 
good quantitative test for constitutive relations of gran- 
ular hydrodynamics. 



IV. MORE RESULTS ON MARGINAL 
STABILITY 

In this Section we investigate the marginal stability 
problem in more detail and obtain some analytic results 
in the dilute limit. 



A. Localization and universality 

Let us characterize the strip state by the scaled param- 
eters 77 and Zq and introduce a different rescaling of the 
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FIG. 4: The critical aspect ratio Ai(/) for the symmetry- 
breaking instability as computed for the constitutive relations 
of Grossman et al. (dash-dotted line) and of JR @ 
(dotted line). Parameter 77 = 11,094. 



coordinate: x — xrj 1 ^ 2 . In terms of the original, physical 
coordinate x p h the new rescaling is independent of the 
system size: 

2n\ 1/2 (l-ry/*x ph 
3 J d 

Equations (g) for the strip state become 

(FZ')' =Q, Z{x = Q) = Z 



and Z'(x = 0) = 0, 
while the eigenvalue problem (|^) and ([?]) reads 



(9) 



= 0, 
and ( 



b'( x = 0) = 



(10) 



Now the primes denote the derivatives with respect to 
x, while k = kr\~ x l 2 is the new scaled wave number. 
Like the coordinate x, the new scaled wave number k is 
independent of L x : 



1/2 



kph d 



(l-r)Va ' 



where k p h is the physical wave number. The problem 
(^) and ( pXj| ) is determined by two parameters: Zq and 
77. However, 77 enters the rescaled equations only in one 
place: in the last boundary condition in Eq. JIo| ) where 
it determines the scaled system size. If the wave func- 
tion 4>(x) is strongly localized in the potential well U(x) 
(correspondingly, the negative energy level is sufficiently 
deep), the results for k and <p(x) become independent of 77 
at sufficiently large 77. Indeed, in this case one can safely 
move the boundary x — 77 1 / 2 to infinity. It is important 
that, in this case, the exact form of the boundary con- 
dition at the driving wall becomes insignificant, leading 
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only to exponentially small corrections [ pl| . This univer- 
sal "localization regime" was discovered in Ref. |l7j that 
employed a different set of constitutive relations. The cri- 
terion for localization can be obtained from the require- 
ment that the localization length (which is of order fc _1 ) 
be much smaller than the (scaled) system size 77 1 / 2 . In 
the physical units it corresponds to the short-wavelength 
limit of the bifurcating solution: k p hL x 3> 1. To fulfill 
this criterion, k should be far enough from the borders 
of the instability interval, where k vanishes. In the next 
subsection we will work in the dilute limit and rewrite 
this criterion in terms of the governing parameters of the 
problem. 

Figure [| shows the marginal stability curves k — k(Zo) 
for different values of 77, obtained numerically. Instead of 
the borders /1 (77) and fa (77) of the instability interval in 
terms of parameter /, the respective borders in terms of 
parameter Z$ appear. One can see that, for large values 
of 77, the marginal stability curves coincide in a wide re- 
gion of Zq not too close to the borders of the instability 
interval. This region corresponds to strong localization. 
Figures ^-|| show the form of the potential (||) and the 
negative energy level — k 2 in three characteristic cases 
(in these figures we returned to the rescaling of the co- 
ordinates and wavenumber by the system size L x ). Fig- 
ure g corresponds to the region of parameters where the 
energy level is deep and eigenfunction is localized. Fig- 
ures ^ and U correspond to the parameter regions close 
to the low- and high-density borders of the instability, 
respectively. There is no localization here. Notice the 
qualitative change in the form of the potential near the 
high-density stability border. Figure |9| shows the respec- 
tive eigenfunctions in these three cases. 



Marginal stability borders and physics of the 
instability 



The low- and high-density stability borders fi(rf) and 
72(77) (or respective borders in terms of Z ) are deter- 
mined by the zero-eigenvalue (k = 0) solution of Eq. (10). 
This solution can be found if we know the strip solution 
Z(x;Z ) of Eq. (g). Indeed, it is easy to check that func- 
tion 0o (x) = FdZ/dZo is a solution of Eq. Jic| ) with 
k = 0, satisfying the boundary condition <f>'(x = 0). Em- 
ploying the second boundary condition cf>(x — 77 1 / 2 ) = 
we obtain dZ 1 /dZ Q = 0, where Z x = Z(x = v 1/2 ;Z Q ). 
For a given 77, this equation is an algebraic equation for 
Zq. Our numerical results imply that this equation has 
only two solutions corresponding to the low- and high- 
density instability borders. The instability borders have 
a clear physical meaning which sheds light on the physics 
of the instability. Let us consider the granular pressure 
p = nT(l + 2G) pfl of the strip state, and introduce a 
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FIG. 5: Marginal stability curves for different values of 77, 
plotted in coordinates k = kr\~ x ^ 2 versus Zo. For a fixed 77 the 
strip state is stable above the respective curve and unstable 
below the curve. The inset shows the splitting of the curves 
near the high-density stability border. The values of 77 are: 
10 4 (solid line), 2.5 • 10 4 (dashed line), 5 • 10 4 (dotted line) 
and 10 (dash-dot line). 




FIG. 6: An example of the negative energy level E — —k 2 ~ 
—83.01 (dashed line) in the potential U(x) (solid line) in the 



regime of localization. The parameters are r/ 
Z = 6. 



2.5 • 10" and 



scaled pressure 
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n c T 



1 + 2G T 

z n 



As P is independent of the coordinates, we can compute 
it at the thermal wall x = 77 1 / 2 . Here T = T and Z = Z\, 
so we arrive at 
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FIG. 7: An example of the negative energy level E = —k 2 ~ 
—0.155 (dashed line) in the potential U(x) (solid line) in the 
absence of localization. The parameters n = 2.5 ■ 10 4 and 
Zo — 75 correspond to the region close to the low-density 
stability border. 
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FIG. 8: Another example of the negative energy level E — 
—k 2 ~ —4.1 (dashed line) in the potential U(x) (solid line) in 
the absence of localization. The parameters n — 2.5 ■ 10 4 and 
Zo — 2.3 correspond to the region close to the high-density 
stability border. 



Now let us compute the derivative dP/df at a constant 
■q: 



dP 
df 



dP dZi dZ Q 



One can easily check that the first and third multipliers 
in the right hand side of this relation are always negative. 
Therefore, the sign of dP/df is determined by the sign 
of dZi/dZo. As we have seen, the marginal stability bor- 
ders are determined by equation 8Zi/dZ = 0. There- 
fore, the steady-state pressure has its extremum points 
points exactly at the points /i and f%. Figure shows 
an example of the dependence P — P(f) at a constant r\ 
for the strip state, found numerically. One can see, that 
P is a decreasing function of / within the instability in- 
terval (/ij/b), and an increasing function of / outside 



FIG. 9: Eigenfunctions <j>(x) corresponding to the eigenvalues 
shown in Figs. ^ (strong localization, solid line), [j] (dash-dot 
line) and ^| (dashed line). The eigenfunctions are (arbitrarily) 
normalized so that 0(0) = 1. 



the interval. The physical interpretation of these results 
is clear. The presence of the anomalous (falling) part 
of the P(f) curve indicates instability, and it is caused 
by the destabilizing role of collisional heat losses @, p"2| . 
We can say that, on the interval (/i,/2), the granulate 
has negative lateral compressibility. At / < /i the heat 
losses are too small to cause instability. The presence of 
the high-density border fa is caused by the finite-density 
corrections to the constitutive relations (that is, by the 
finite size of the particles). This is in contrast to radia- 
tive condensations in gases and plasmas, where such a 
stabilizing effect would be absent. Now consider a small 
density modulation of the strip state with a very long 
wavelength: k — > 0. For this perturbation, the stabiliz- 
ing effect of the lateral heat conduction vanishes, and the 
negative compressibility makes the strip state unstable on 
the interval (/i,/2). For a non-zero k, the lateral heat 
conduction has a stabilizing effect. Therefore, a density 
modulation with too a short lateral wavelength should 
be stable, as Fig. || indeed shows. 



C. Dilute limit 

In this subsection we shall work in the dilute limit and 
derive several analytic results. We shall see that, at large 
rj, the dilute limit faithfully reproduces the low-density 
parts of the marginal stability curves. 



1. Strip state and marginal stability problem 

In the dilute limit, Z> 1, Eq. (0) for the strip state 
becomes 



(Z 3/2J// = 3Z -l/2 



(11) 
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FIG. 10: The scaled steady-state granular pressure P versus 
the grain area fraction / for the strip state. The two circles 
correspond to the marginal stability borders /i and fa. The 
effective lateral compressibility of the gas is negative on the 
interval (/i,/2). Parameter rj = 10 4 . 



where the primes denote the x derivatives. The boundary 
conditions are Z'{x = 0) = and Z(x — 0) = Zq, where 
Zq is related to / and rj by the normalization condition 
fff* Z~ 1 dx = jr] 1 / 2 . The solution of this problem is 
elementary: 



= -y (arccoshV? + V C 2 - C 



where £ = Z/Zq and 



Zn 



27? 



1/2 



f v i/2 + i s i n h(2/r; 1 /2) 



(12) 



(13) 



(Returning for a moment to the old rescaling of the coor- 
dinate, x p h/L x — > x, one can see that the density profile 
( |l2| ) is determined by a single parameter: £ = ft] 1 / 2 .) 
Eq. (|I^) for cj>{x) = (%/3/2) Z x t 2 {x) ip k (x) takes the form 



1 

c 2 " 



0. 



(14) 



where k = k 2 Z$, and £ = C(x) is given, in an implicit 
form, by Eq. (p"2|). The boundary conditions for Eq. ( p^ ) 
remain the same as in Eq. (]To|) . As C(^) is a monotonic 
function of a;, we can change the independent variable in 
Eq. ( ful ) from 5 to £. The resulting equation for <ft(() is 



4(C-l)C0" + 2^' + (l~«C 2 )^ = O, 



(15) 



where the primes now denote the £ derivatives. The 
function </>(() is defined on the interval 1 < C < £i, 
where Ci = Zi/Zq = cosh 2 £. One boundary condition 
is </>(£ = £i) = to which we may add an arbitrary 
normalization condition (/>(£ = 1) = 1. An additional 
boundary condition, <fi' {Q = 1) = [k — l)/2, can be found 
from Eq. (Il5|) itself, after substituting there Q = 1. This 



eigenvalue problem includes a single parameter £, while 
k serves as the eigenvalue. 

We have been unable to solve Eq. ( |l5| ) analytically for 
a general k. Still, several important asymptotics can be 
obtained. 



2. Zero-energy state and stability border /i 
For k = Eq. (|l|) can be solved analytically: 

0(C, k = 0) = MO = VC ~ VC~l arccoshv/C, (16) 

In other words, we impose a zero eigenvalue k = and 
find the low-density stability border fi = fi (rf) from the 
boundary condition </>o(Ci) = 0. We obtain an algebraic 
equation coth(£i) = £i for £i = fiij 1 / 2 . Its solution is 
£i = 1.19968 .... This result explains why all marginal 
stability curves shown in Fig. [| depart (almost) from the 
same point at the low-density side. Figure |ll| compares 
the scaled quantity fir] 1 / 2 at different rj, found numeri- 
cally from Eq. (^), with this analytic prediction (a con- 
stant). The agreement is very good for large rj. As rj goes 
down, fi increases and the dilute approximation starts 
to deteriorate. 




FIG. 11: The scaled low-density stability boundary £i = 
fin 1 ^ 2 at different rj 1 ^ 2 as found numerically from Eq. (^) (cir- 
cles). Solid line shows the analytical result f i = 1.19968 . . . 
obtained in the dilute limit. 



Figure [12] compares the analytic result for the zero- 
energy eigenfunction 4>q(x), given by Eqs. ( |l6| ) and (|l2j), 
with a numerical solution of Eq. (|6|) for rj — 10 4 . The co- 
ordinate x in Fig. [l2| is rescaled by L x . The analytic and 
numeric results are obtained for slightly different values 
of /i (see Fig. |ll]). One can see that the agreement is 
excellent. 
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FIG. 12: The zero-energy eigenfunction (j>o(x) computed nu- 
merically from Eq. (ph (circles) and given analytically by 
Eqs.(^) and ( |l^ ) (solid line). The parameters are r\ = 10 4 , 
fi,num = 0.0124 (circles) and fi,anai = 0.0120 (solid line). 



3. Short-wavelength limit: localization and universality 

In the short-wavelength limit the system boundary 
C = Ci = cosh 2 £ can be moved to infinity. This re- 
quires a strong inequality £ = frj 1 ' 2 1. In this limit, 
the eigenvalue problem ( |i~5j ) does not include any param- 
eter. The eigenvalue k should therefore be a number of 
order of unity, hence k = A/Z , with constant A of or- 
der of unity. The constant can be found numerically: 
A ~ 0.525. This simple result represents the low-density 
limit of the "universal" marginal stability curve, corre- 
sponding to strong localization. Figure |l3| shows this 
asymptotics for r\ = 10 5 . One can see excellent agreement 
for large enough Zq, but not too close to the higher- Zq 
(low-density) instability border. Near the instability bor- 
der £ becomes of order of unity, and localization breaks 
down. 

0.06F 



0.04 



s 



0.02 




FIG. 13: The marginal stability curve for n = 10 5 , plotted 
in coordinates k = krj" 1 ^ 2 versus Zo (solid line) and two 
dilute- limit asymptotics: the dilute- limit part k = A/Zo of 
the ( "universal" ) short-wavelength curve (dotted line) and the 
long- wavelength asymptotics fl24) (dashed line). 



Returning to the parameters rj and /, and to the wave 
number k = k p hL x , we can rewrite the asymptotics k = 
A/Z as 



.4 



iV/2 + \ sinh^/r? 1 / 2 ) 



(17) 



This asymptotics is shown in Fig. FT 
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FIG. 14: Marginal stability curves for different values of 
?7, plotted in coordinates k versus i = fv 1/2 - For a fixed 
rj, the strip state is stable above the respective curve and 
unstable below the curve. The values of n are: 10 4 (solid 
line), 2.5 ■ 10 4 (dashed line), 5 ■ 10 4 (dotted line) and 10 5 
(dash-dot line). Also shown is the dilute-limit asymptotics 
( p^| ) of the universal marginal stability curve. 

The asymptotics ( fl7| ) is valid when £ = ^> 1. 

It is easy to check that this criterion coincides, in the 
dilute limit, with the localization criterion k^ 1 ry 1 / 2 
discussed in Sec. 4A. On the other hand, the parameter £ 
should not be too large, so that the dilute limit condition 
Z Q 1 is still satisfied, see Eq. (p"3|). These two criteria 
can be rewritten as a strong double inequality for /: 

1 . log^T? 1 / 2 ) 

< / < 



1/2 



n 



2^/2 



that can be satisfied only for extremely large rj. 



4- Long-wavelength limit: perturbation theory 

Close to the low-density stability border, / — /i <C 
fi = 1.19968 .. .v^ 1 !" 1 we can assume that k -C 1 and 
solve Eq. (|l5|) perturbatively. In the physical units, 
this strong inequality corresponds to the long-wavelength 
limit: k p hL x <C 1. In its turn, the dilute limit requires 
Tf > 1. We substitute in Eq. (ph 0(C) = MO + «^(C), 
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where 4>q is given by Eq. ( |l6| ) . Neglecting the K 2 -term in 
Eq. (|l5|), we obtain 



4(c - 1) C^" + 2^' + n = c 2 <MC) 



(18) 



The normalization and boundary conditions at £ = 1 are 
ft(C = 1) = and fi'(C = 1) = 1/2, respectively. The 
latter condition follows from Eq. ( p!s|) itself. Equation 
( |l8| ) can be solved analytically. With the account of the 
two boundary conditions we obtain 



i 



M0 + M0h (C)-*i (0^2 (0, (i9) 

-20o (0- /i and 



where $i(C) = (C - 1) 1/2 and $ 2 (C) 
^2 are indefinite integrals: 



h 



£iGi 



g?£ and I2 



$2Gi 

G 2 W 



dC, (20) 



where d = C 2 0o/4 , G 2 = C(C - 1), W = <J>i$ 2 - $i$2 
and the primes denote the C derivatives. Integrals I\ 
and I2 can be evaluated in elementary functions, but the 
results are too cumbersome to be presented here. The 
additional boundary condition [at ( = Zi/Zq = £1] reads 
= ^o(Ci) + «0(Ci) = which yields 



«(0 = -0o(Ci)/n(Ci) • 



(21) 



where £1 = cosh 2 £. At the low-density stability border 
/ = /1 we have £ = £1 = 1.19968 .... In the perturbative 
treatment, one should expand near £ = £1 up to the 
linear term £ — £1. The zero-order term vanishes, and we 
obtain 



«(0 



2$i(£i)I a (£i) 
In the physical variables we have 

~ 2$ 



i(a)/ 2 (a)y ^ ?u 



(22) 

As Zo is proportional to 77 1 ' 2 , the right hand side of Eq. 
( p2| ) is actually independent of 77. Evaluating the integral 
we obtain 7 2 (£i) = -0.883381 . . .. The final result 

is 



kphL x = 2.5115 
Alternatively, we obtain 

= 0.000485 



1/2 



(23) 



3/2 



Zp - gbfo) 
2/2 (£1) 
(Zo^O-Zo) 1 ^ 



1/2 



(24) 



The asymptotics (24) is depicted in Fig. 03. Close to the 



higher-Zo (low-density) stability border it shows good 
agreement with the marginal stability curve found nu- 
merically. 



V. SUMMARY AND DISCUSSION 

We determined the criteria for the spontaneous 
symmetry-breaking instability of the laterally uniform 
granular cluster (strip state) in a prototypical driven 
granular gas. Working in the limit of nearly elastic parti- 
cle collisions and low or moderate densities, we employed 
granular hydrodynamics with the Jenkins-Richman con- 
stitutive relations p0| . The instability of the strip state 
can be interpreted in terms of negative compressibility 
of the granulate in the lateral direction. An important 
limit is found, when the marginal stability curves are in- 
dependent of the details of the boundary condition at 
the driving wall. In this regime the density perturba- 
tion is exponentially localized at the elastic wall opposite 
to the driving wall. Working in the dilute limit, we ob- 
tained some analytic asymptotics of the marginal stabil- 
ity curves. 

The results of this work show that the symmetry- 
breaking instability predicted in Ref. jL7| is robust and 
does not require very special constitutive relations. The 
marginal stability curves obtained in this work are quite 
similar to those obtained earlier |l7j for a different set of 
constitutive relations (see Fig. ||). There are some quan- 
titative differences, however. Therefore, the instability 
provides a sensitive test to the accuracy of constitutive 
relations. 

This work was focused on the criteria of instability 
of the strip state. In systems sufficiently long in the 
lateral direction, instability occurs in a whole range of 
wave numbers k (below the respective marginal stability 
curve). Correspondingly, multiple steady state solutions 
with different k are possible. In a laterally infinite sys- 
tem, these solutions are periodic in the lateral coordinate. 
A finite system selects a finite number of wavelengths 
p7{ . An important issue that was not addressed in this 
work is selection: what is the wavelength of the resulting 
symmetry-broken cluster in an infinite, or long enough, 
system? The selection has dynamical nature; this impor- 
tant issue is addressed elsewhere p2f . 

Recently, the predicted symmetry-breaking instability 
has been observed in particle simulations psj . We hope 
it will be investigated in experiment, too. The experi- 
mental setting can be of the type used by Kudrolli et al. 
0, [r|: a system of steel spheres, rolling on a smooth 
surface and driven by a rapidly vibrating side wall. The 
present work (see also Ref. ]l^] ) provides the region of 
parameters where the instability can be observed. An im- 
portant issue is to eliminate the static friction between 
the particles and surface that occurs far enough from the 
driving wall. In experiment, this is achieved by slightly 
inclining the system, so that a very small gravity appears 
j7, 16 . As the result, the strip state moves down, toward 
the driving wall 0. The model problem investigated in 
the present work does not include gravity. We expect, 
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however, that the symmetry-breaking instability should 
persist for a non-zero gravity. In fact, a similar insta- 
bility has already been observed in particle simulations 
of a dilute two-dimensional granular bed fluidized by a 
rapidly vibrating bottom plate 1 24 . Under conditions of 
the simulations [S3 there was no direct, mechanical cou- 
pling between the bottom plate vibration and collective 
granular motions. Therefore, the vibrofluidized system, 
investigated in Ref . |24| , is similar (though not identical) 
to the model system driven by a thermal wall. As grav- 
ity introduces an additional scaled parameter, the phase 
diagram of this type of systems should be more compli- 
cated. For example, it is already known that, at some 
values of the scaled parameters, steady "thermal" con- 
vection (steady state of a different type) develops both 
in vibrofluidized systems J24|, and in systems driven 
by a "thermal" wall |26|, |27jj. Granular hydrodynamics 
will be instrumental in delineating the phase diagrams of 
these systems in the limit of nearly elastic collisions. 

Finally, when inelasticity of the particle collisions is 
not small, the normal stress difference, possible lack of 
scale separation and non-Gaussianity in the velocity dis- 
tribution may become important. The potential role of 
these effects in the symmetry-breaking instability should 
be the subject of further investigations. 
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